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Abstract 



CN ' We construct the explicit two-holon eigenstates of the SU(2) Kuramoto- Yokoyama model at 

the level of explicit wave functions. We derive the exact energies and obtain the individual holon 
^ , momenta, which are quantized according to the half-Fermi statistics of the holons. 
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I. INTRODUCTION 



One milestone towards the understanding of fractional quantization in one dimension 
is the model independently introduced by Haldane Q| and Shastry P| in 1988. The 
model describes a spin 1/2 chain with a Heisenberg interaction which falls off as one over 
the square of the distance between the sites. The exact ground state is provided by a 
trial wave function proposed by Gutzwiller as early as in 1963. The Haldane-Shastry 
Model (HSM) offers the opportunity of studying spinous, i.e., the elementary excitations 
of one-dimensional spin chains, on the level of explicit and analytical expressions for one 
and two-spinon wave functions 0|, which are at least at present not available for any other 
model. Kuramoto and Yokoyama generalized the model to allow for mobile holes [i.e., 
empty lattice sites) with a hopping parameter that also falls off with 1/r^ as a function 
of the distance. The Kuramoto- Yokoyama Model (KYM) hence contains spin and charge 
degrees of freedom, and accordingly supports spinon and holon excitations, which carry spin 
I but no charge and charge +1 but no spin, respectively. In principle, the KYM allows for a 
similarly explicit construction of holon wave functions, which so far have only been obtained 
for states involving a single holon. The reason for this deficit has been of technical nature, 
related to the commutation relations of the operators used to build the Hilbert space of 
these fractionally quantized excitations. Whereas for the spinous one can use bosonic spin- 
flip operators, one needs fermionic creation and annihilation operators for the holons. 

In this article we address and overcome this technical problem as we construct the explicit 
wave functions for two-holon excitations of the KYM. The article is organized as follows: In 
Section ITTl we review the KYM and its properties. In Section ITTTl and HV| we briefly discuss 
the ground state at half filling and the spinon excitations. We further review the analytic 
results so far known for the one-holon excitations in Section |3 as a preliminary for the 
construction of the explicit two-holon wave functions to be done in Section IVIl Therein we 
derive the exact energies and individual holon momenta, which turn out to be quantized 
according to half-Fermi statistics of the holons. 
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II. KURAMOTO YOKOYAMA MODEL 

The Kuramoto-Yokoyama model P| is most conveniently formulated by embedding the 
one-dimensional chain with periodic boundary conditions into the complex plane by mapping 
it onto the unit circle with the sites located at complex positions 77^ = exp(i^a), where N 
denotes the number of sites and a = 1, . . . , N. The sites can be either singly occupied by 
an up or down-spin electron or empty. The Hamiltonian is given by 

where Pajj exchanges the configurations on the sites rja and 77/3 including a minus sign if 
both are fermionic. Rewriting (0) in terms of spin and electron creation and annihilation 
operators yields 



H - V ^ P 
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Pg, (2) 



where the Gutzwiller projector Pq = Y[ai^'~'^a^'^ai'^ai'^a^) enforces at most single occupancy 
on all sites. The charge occupation and spin operators are given by Ua = cl^iC^-^ + cl^^c^^ 
and = ^ J^aa' '^aa'^'^cr' ^ao-' ^ whcrc t", = x,y, z, denote the Pauli matrices. 

The interaction strength in is an analytic function of the lattice sites by use of 

1 _ VaV, (3) 



The KYM is supersymmetric, i.e., the Hamiltonian commutes with the operators 
a^aa^ab^ whcrc ttaa dcuotcs the annihilation operator of a particle of species a (a 
runs over up- and down-spin as well as empty site) at site rja. The traceless parts of the 
operators J"'' generate the Lie superalgebra su(l|2), which includes in particular the total 
spm S = J2a=i Sa- In addition, the KYM possesses a super- Yangian symmetry [G], which 
causes its amenability to rather explicit solution. 



III. VACUUM STATE 



We first review the ground state at half filling, which is the state containing no excita- 
tions (neither spinous nor holons). For even, this vacuum state is constructed by the 

3 



Gutzwiller projection of a filled band (or Slater determinant (SD) state) containing a total 
of electrons: 

\^o)=Pg n iAl\^)^Po\^^u)- (4) 

|9|<'?F 

Taking the fully polarized state |0|) = HqcIjiIO) as reference state, we can rewrite the 
vacuum state as 

\^o) = 5^ ^o(^i, . . . , zm) St, . . . S+„ |0;) , (5) 

where M = N/2 and the Zi^s denote the up-spin coordinates. The sum in (jSJ extends over all 
possible ways to distribute the coordinates Zi's over the lattice sites The wave function 
is given by 1^ 

M M 

"^oizi, ...,zm) = Ylizi - Zjf Yl (6) 

i<j i=l 

its energy is 

= - — . (7) 

The total momentum of a state is evaluated by considering the operator T, which translates 
all coordinates counterclockwise by one site. T is related to the momentum operator P via 

T = exp(-zP). (8) 

This yields the momentum of |\l'o) to equal zero if is divisible by four and vr otherwise. 

Note that (jHI) represents the ground state of only at half filling, i.e., when all sites 
are occupied. As was shown by Kuramoto and Yokoyama j^, the ground state away from 
half-filling can be constructed by Gutzwiller projection similar to (jH). 



IV. SPINON EXCITATIONS 

Let A^ be odd and M = {N — l)/2. A localized spinon at site "r^^" is constructed by the 
Gutzwiller projection of an electron inserted in a Slater determinant state of A^ + 1 electrons: 

|*^''> = ^GC,,K+^>. (9) 

The annihilation of the electron causes an inhomogeneity in the spin and charge degree of 
freedom. After the projection, however, only the inhomogeneity in the spin survives. The 
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spinon hence possesses spin one-half but no charge. The wave function of a locahzed spinon 
is given by B 

M 

^^P(2i, ...,zm) = n(^7 - ^0 ^0(^1, ■■■.Zm), (10) 

i=l 

where \l/o is defined in (jSJ. Fourier transformation yields the momentum eigenstates 

N 

i« = E(^7r|^T>' (11) 

a=l 

which vanish identically unless < m < M. In particular, this implies that the localized 
one-spinon states form an overcomplete set. It is hence not possible to interpret the 
"coordinate" t]^ literally as the position of the spinon. The momentum eigenstates (fTTjl are 
found to be exact energy eigenstates of the KYM, with its energies given by B 

E'^ A (12) 



m \ 2 

The spinous obey half-Fermi statistics, which was first found by the investigation of their 



spinouf 

fin 



state counting rules jS]- Later it became apparent that the fractional statistics of the spinous 
manifests itself in the quantization rules for the individual spinon momenta as well 

V. ONE-HOLON EXCITATIONS 



The charged elementary excitations of the model are holons, the concept of which must 
be invoked whenever holes and thereby charge carries are doped into the chain. A localized 
holon at lattice site 77^ is constructed as 

l^r>=%^G4Jv^^D^>. (13) 

(Alternatively we could use the operators c^^ and c^^.) Compared to the spinon we elim- 
inate the inhomogeneity in spin while creating an inhomogeneity in the charge distribution 
after Gutzwiller projection. Thus the holon has no spin but charge e > (as the electron 
charge at site r]^ is removed). Note that the holon is strictly localized at the holon coor- 
dinate ^, as holon states on neighboring coordinates are orthogonal. In total, there are N 
independent one-holon states (fTT^ . 



Momentum eigenstates are constructed from p3|) by_ 



that only [N + 3)/2 of them are energy eigenstates 



burier transformation. It turns out 
. We will restrict ourselves to this 



subset in the following. These states are readily described in terms of their wave functions. 
We take |0|) as reference state, and write the one-holon energy eigenstates as ^| 

K°> = E (^1' ■■■^^M;h) c,^ . . . Si, |0i) , (14) 

{Zi;h} 

where the sum extends over all possible ways to distribute the up-spin coordinates Zi and the 
holon coordinate h over the lattice sites rja subject to the restriction Zi 7^ h. The one-holon 
wave function is given by 

M 

. . . , ZM\ h) = - z,) ^0(^1, . . . , ^m), (15) 

i=l 

where is given by (jH))- Note that as a sum over the coordinates h is included in (I14|) . no 
such sum is required in (jl5j) . It can be shown that < m < M + 1, where M = {N — l)/2 
is the number of up-spin coordinates, the wave function (fT^ represents an exact energy 
eigenstate with energy [ill 

For other values of m, the states j^E'J^) do not vanish identically (as I^E'^) for spinous do), 
but are not eigenstates of the Kuramoto-Yokoyama Hamiltonian either. Consequently, 
we are allowed to refer to the states (fT^ with (fT3j) as "ho Ions" only ifO<m<M + l. 

Note that this implies that the states (fT^ do not really constitute "holons" localized in 
position space, but only basis states which can be used to construct holons if the momentum 
is chosen adequately. The total number of single-holon states is given by M-l-2, according to 
the number of distinct values m is allowed to assume. Since the states (jl3|) are orthogonal 
for different lattice positions ^, there are N = 2M + 1 orthogonal position basis states 
|\E'^°). Hence the states |^E'^°) cannot strictly be holons, but rather constitute incoherent 
superpositions of holons and other states. It is clear from these considerations that it is not 
possible to localize a holon onto a single lattice site. The best we can do is to take a Fourier 
transform of the exact eigenstates j^^J^) for < m < M + 1 back into position space. The 
resulting "localized" holon states will be true holons but will not be localized strictly onto 
lattice sites. Such a true holon state "localized" at a given lattice site will not be orthogonal 
to such a state "localized" at the neighboring lattice site, as there are only M + 2 holon 
states while there are lattice sites. The situation is hence very similar to the case of 
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the spinons, which form an overcomplete set and are well known to be non-orthogonal if 
"localized" on neighboring lattice sites. 

The one-holon momenta of the states (|14j) with (|15|) are derived in analogy to the vacuum 
state to be 

p'^ = In+^[^~\) mod2vr. (17) 
If we introduce the one-holon dispersion 

we can rewrite ()16|) with the vacuum energy ((7j) as 

E^ = Eo + e^°(p::,°). (19) 

VI. TWO-HOLON EXCITATIONS 
A. Momentum eigenstates 

Let be even and M = {N — 2)/2. The two-holon state with holons localized at r^^^ and 
T]^^ is constructed in analogy to (fT3|) as 

\^t,.)=%i%iPocl,clJ^^^-'). (20) 

In analogy to ()15|1 . a momentum basis for the two-holon eigenstates is provided by the 
wave functions 

M 

<°„(^1, . . . , Zm; h, h2) = {h - /l2)(/ir^2 + ^l^D 11(^1 - - Z,)^o{z,, ...,Zm), (21) 

1=1 

where \&o is again given by (jHI), /ii,2 denote the holon coordinates, and the integers m and n 
satisfy 

0<n<m<M + l. (22) 
The corresponding state is then given by 

{zi;hi,h2} 
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where the sum extends over all possible ways to distribute the up-spin coordinates Zi and 
the holon coordinates hi^2 over the lattice sites rja subject to the restriction Zi ^ h\ ^ h2. 
The momentum of the states is easily found to be 

p!^°n = ^Ar+^(m + ri) mod27r. (24) 

It can further be shown that the states (|23|) are spin singlets, i.e.^ they are annihilated by 
as well as . 



B. Action of -f^KY on the momentum eigenstates 

In the following we will construct the two- holon energy eigenstates starting from ()2ip. 
First, we introduce the auxiliary wave functions 

M 

^mn{Zi, ...,Zm; hi, h2) = K^h^ W{hi - Zi){h2 - Zi) ^q{zi, Zm)- (25) 

2=1 

The action of the Hamiltonian on the states (j2T| will be obtained later via 

^mn — ^m+l,n + V^n.+l,m ~ V'm.n+l ~ V^n,m+1- (26) 



Second, we rewrite the Hamiltonian in analogy to as 

Hky = ^ [hi- + H'i + Hy + Hl + H^) , (27) 

where we separate the spin-exchange, spin-Ising, potential, j-charge kinetic term, and |- 
charge kinetic terms. In the following we treat each term separately. 

For the spin-exchange term we begin by observing that [S^Sp(pnm]{zi, . . . , zm', hi, /i2) is 
identically zero unless one of the arguments zi, zm equals rja- We have 



H^ifmn {zi,. . . , Zm', hi, h2) = 



/ 1 12 '■Pmn 



{zi, . . .,zm; hi, h2) 



Sr^Sr^ Vl3 ^nm[Zl, ■ ■ ■ , Zj-i,r]f3, Zj+i, . . . , Zm'-i III, 112) 

i=i l^^j I ^ ''^i 

f Vf ^/3(r/,-..)' ( ^U^i,...,ZM;hi,h2) 

\ ^ \ ^ , O (pmn 

j=i 1=0 J ■' 



M 



^ , (iV-l)(iV-5) ^ _ N-3 ^ 1 \ if 



12 



<^ — (5 - 2N)hTh^, - j2 1 _ 12 - ^r^r' 



fj2 f)2 
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2 dh, 



3 \ hmun+2 o 7,^+2/1" 1 ,o 



/ hi — h2 dh2 hi — /i2 dhi [ /i^/ij 



(28) 



where we have introduced the coefficients Ai = — Yla=i VaiVa — 1)' ^- Evaluation of the 
latter yields Aq = (N - 1){N - 5)/12, Ai = -{N - 3)/2, A2 = 1, and A, = for 2 < / < 
A^- 1 Q. 

For the spin-Ising term we obtain 

N 



(zi, . . . , Zm] hi, h2) 
1 



M M 

2 l^j — hir ^ \zi — h'} 



M 



4 = 1 



(^1, • • • , 2;^,/ ; ^1, ^2) 

|/ll-/l2P 



_ AT- 



48 



The potential term yields 



(29) 



Hy^mn (^l5 • • • 5 ^Af! ^2) = 
1 1 



E-Pg {-\nonfi + - I) -Pg 



_ 1 A^2 _ ^ 



(zi, . . . , zm; ^2) 



(30) 



2|/ii-/i2p 12 4 12 

The charge kinetic terms deserve particular care as new techniques are required. 
For the |-charge kinetic term, we ffist observe that [c^|C^|V9m„](2;i, zm\ hi, /i2) is iden- 
tically zero unless one of the arguments zi, ...,zm equals rja- We thus find 



H^iPmn (Zi, . . . , Zm] hi, h2 



■ N 

E 



\Va-V/3? 



{zi, . . . , Zm', hi, h2) 



a=hi,h2 Pf^a 
N 

E 



[Zi, . . . , ZM'ihi,r]p) s. {zi, Zm] ^7/3, h2) 



+ E 
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^ ^ V^^iVp - ^2)' & ( iPmn{Zl, ■■■,Zm; hi,T]p) 



1=0 P^h, ^'1/^2 9r/^ 
M N 



/=0 I3j^hi 



l\\hi - r]p\^ dri^p I r]f 



iV^ - 1 m(m - A^) n{n - A^) 



6 

N -\ 



+ 



9^ 



iV- 1 



n 



_d_ 



I ^y^m+2,n , ^ . mr n+2 1 • • • , hi, /j-2) 



(31) 



where we have introduced the coefficients Bj^ = — Yl'^=i V'fs^^iv^^^Y^^^ which are evaluated 
in Appendix 1X1 (p?T|l is vahd if and only if < n,m < (A^ + 2)/2, which finally leads to the 
restriction (j22I) for the actual ^E'^n's- 

For the treatment of the J,-charge kinetic term we avail ourselves of the fact that ipmn can 
be equally expressed by the up-spin or down-spin variables, as we show in Appendix |B| If 
we denote the down-spin coordinates by Wi, we obtain 



H^fmn {Zl, • • • , Zm', hi, h2 



N 



{zi, . . . , Zm] hi, /i2) 



Ar2 



6 

A^- 1 



+ 



m(m — A^) n{n — A^) 



d_ 

dhi 



dh,. 



, l;,m+2.n , ^ 1 Vmn{wi, . . . , W M] hi, 



dhi - 2^''^' dh. 



(32) 



Using identities verified in AppendixOfor the derivatives with respect to the Zj's and /ii,2's, 
the total charge-kinetic term becomes 



Ar2-1 



(2:1, . . . ,2;m; ^1, h'2) 



+ m(m - A^) + n(n - A^) - C2 - + m Ci - 



1\ hi + h2 fm-n 



92 



92 
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{hi — dh2 ^ h2 — hi dhi 



hi 



-Cihl 

M 



+ 



dh_ 

V-'mn (2^1 



{hi - h^Y 



{Zi - IL2 
Zm', hi, /l2 



hi 



d_ 

dhi 



{zi - hy 



(33) 



with the constants Ci = 1/(1 - = (N - l)/2 and C2 = ^2=i V(l - Va? = 

{6N — 5 — N'^)/12 introduced and evaluated in I4|]. Summing up all terms, we finally obtain 
the action of the Hamiltonian on the auxiliary wave functions ipmn'- 
27r2 



^7V-l 



iV2 

l hi + h2 
2hi — h2 

With this implies 
27r2 



8-9Ar . . N-2 
h m(m — JM) + n(n — A/ j + m h n 



ho 



Ar2 



n) + 



N 



hi 



hi 



{hi - h2y 



(34) 



m m+ n nN/ 



,ho 



m — n 



{hi 



/^2)^^(/^r/^2-/^i/^^)^o 



27r2 



N 



hi — h2 



N 



\- \ m ] m + \ n ] n 



N 



m — n 



ho 



(m — n) \E' 



ho 

m-— 1,71+1'' 



1=1 



where we have used I 



x—y 



rr-m-l„,n+l 



X 



y 



(35) 



(x™y" + x"y™). The symbol [J 



denotes the fioor function, z.e., [xj is the largest integer / < x. First, note that the action 
of the Hamiltonian on the \I/J^„'s is trigonal, i.e., the "scattering" in the last line is only to 
lower values oi m — n. Second, shows that the states form a non-orthogonal set. 
We will now proceed to construct an orthogonal basis of energy eigenf unctions. 



C. Energy eigenstates 

Due to the trigonal structure of the Hamiltonian when acting on the ^l/^n's we can derive 
the energy eigenstates using the Ansatz 

|<^'::^°n>= E «r|^™%n+/>, (36) 
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which yields the recursion relation 

^ a^"(m - n - 2k), a^" = 1. (37) 



1 



2l{l -l + n-m) 



k=0 



This defines the two-holon energy eigenstates (jHBj) . The energies are given by 



where the momentum quantum numbers satisfy 



/ m-n 
m m + \ n nH 



(38) 



< n < m < y , (39) 



and the total momentum is given by 

For the lowest energy state, (jHHj) simplifies (up to an additive constant '7r^/12A^) to the 
ground-state energy of the chain doped with two holes, which is a special case of the result 
by Kuramoto and Yokoyama for the ground state at general filling fraction. 



VII. FRACTIONAL STATISTICS 

Fractional statistics in one-dimensional systems was originally introduced by Haldane 
in the context of non-trivial state counting rules. Recently, it was realized that the fractional 
statistics of spinous in the HSM also manifests itself in specific quantization rules for the 
individual spinon momenta We now apply this line of analysis to the holon excitations 

in the KYM. 

To begin with, let us recall that the asymptotic Bethe ansatz solution of the KYM 0, fl^ 
implies that the holons are free, i.e., that they interact only through their statistics, while 
there is no position or momentum dependent interaction potential between them. This 
induces us to rewrite the two-holon energy (jHHj) as 



E^°„ = Eo + 6'^°(p™) + e^°(p„), (40) 

where we assume the one-holon dispersion (|TH|) and introduce the single-holon momenta 
according to 

-K 2tx ( \\ TT 27r / 1\ 

^- = -2 + ivr+i ' ^" = -2 + ivr-4 • 
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Note that the fractional shifts of ^ ■ | in and p„ occur in opposite directions. Since 
n < m, the momenta are shifted away from each other, implying > Pn- The shifts 
directly follow from ()4()|1 : any other assignment of the single particle momenta would yield 
an additional interaction term in the energy, corresponding to the last term in (|HH|) above. 
For the difference in the individual holon momenta we hence obtain 

Pm-Pn = ^ [- + integer . (42) 



N \2 

We interpret this result as a direct manifestation of the half-Fermi statistics of the holons, as 
the shift in the single particle momenta can be attributed to a statistical phase acquired by 
the states as the holons pass through each other 1^. Indications of the half-Fermi statistics 
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li of the 



of the holons have previously been observed in thermodynamical quantities 
KYM as well as the electron addition spectrum ^5 Q| • 

Let us now elaborate on the general implictions of this result. The wave functions we 
have obtained above are of course eigenstates of the KYM only, which is as idealized as 
integrable and exactly soluble models tend to be. The quantization rules for the single 
particle momenta we have obtained for this model, however, have a much broader validity. 
As mentioned above, the unique feature of the KYM is that the holons are free in the sense 
that they only interact through their fractional statistics. The single particle momenta of 
the holons are hence good quantum numbers, which assume fractionally spaced values. For 
two holons, these are given by ()41|) . The crucial observation in this context is that the 
statistics of the holons is a quantum invariant and as such independent of the details of the 
model. This implies directly that the fractional spacings are of universal validity as well. If 
we were to supplement the model we have studied by a potential interaction between the 
holons, say a Coulomb potential, this interaction would introduce scattering matrix elements 
between the exact eigenstates we obtained and labeled according to their fractionally spaced 
single particle momenta. These momenta would hence no longer constitute good quantum 
numbers. The new eigenstates would be superpositions of states with different single particle 
momenta, which individually, however, would still possess the fractionally shifted values. In 
other words, looking at the quantization condition P2|l . the "1/2" on the left of the equation 
will still be a good quantum number, while the "integer" will turn into a "superposition of 
integers" in the presence of an interaction between the holons. 
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VIII. CONCLUSIONS 



In this article we have studied the two-holon states of the Kuramoto-Yokoyama model. 
We constructed the explicit two-holon wave functions and derived their momenta and ener- 
gies. The results display the half-Fermi statistics of the holons, which manifests itself in a 
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APPENDIX A: B-SERIES 

Evaluation of the series 



shift of 1^ in the difference of the individual holon momenta. 



N-l 




(Al) 



with / restricted to < / < M = (iV 



2)/2 yields 




for < n < A^, 



(A2) 




for < n< A^, 



for n — N, 



(A3) 




1 for < n < A^ - 2, n = AT, 

1-N for n = AT - 1, 



(A4) 




for / > 3, < n < 



A^ + 2 
2 



N + 2 
2 



(A5) 



for / > 3, 



<n< N. 
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Proof: Bq, B^, and B2 are found by straight forward evaluation of the respective sums. For 
i5l) consider 

N-l 1-2 

k 

N 



a=l k=0 
1-2 



k=0 ^ ^ a=l 



k=0 
, n 



5^ \ k 

k=0 ^ 

u—n \ / 



\n—k 



for 3 < /,0 < n < (A^ + 2)/2, 



N5k,N-i~n = N[ / ^1 otherwise. 
A/ — 72 — 1 . 



V fc=0 

For the last steps note that the the binomial coefficients of even and odd sites equal each 
other. 



APPENDIX B: WAVE FUNCTION IN j- AND i-SPIN COORDINATES 



The wave functions '^^^ can be equally expressed either in up-spin (z) or down-spin 
coordinates [w): 



K^K^i^hi - hi) 



M{M+1) Uj (hi - Zj){h2 - Zj)Zj Ui^^i^i - ^j) Ul,jiwi - ^l) 



^mn{Wu-,WM;hi,h2) 



(Bl) 



h'^h^ih - hi) 

This identity applies to the auxiliary wave functions ifmn-, as the prepolynomial contains 
only the coordinates h\^2- 



APPENDIX C: A DERIVATIVE IDENTITY 



The necessary relation for the J,-charge kinetic term is 



E 



hi 



{Wi - h2){Wj - /12) 



hi 



^ (zi - /lor 



E 



h2 h2 



+ 2 



h\ 



ih - h2y 



15 



with Ci and C2 defined as above. ()C1|) is also valid for hi /i2- 
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